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Abstract. 

Specific activator and repressor transcription factors which bind to 
specific regulator DNA sequences, play an important role in gene ac- 
tivity control. Interactions between genes coding such transcripion 
factors should explain the different stable or sometimes oscillatory 
gene activities characteristic for different tissues. In this paper, the 
dynamic P53-Mdm2 interaction model with distributed delays and 
weak kernel, is investigated. Choosing the delay or the kernel's coeffi- 
cient as a bifurcation parameter, we study the direction and stability 
of the bifurcating periodic solutions. Some numerical examples are 
finally given for justifying the theoretical results. 
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1. Introduction 



P53 is a very important gene in oncogenesis. It is also known as " Guardian 
of the genome". Its anomalies are almost universal in tumoral cells [5]. 
The full activity of p53 gene starts when is detected a DNA damage [3,4]. 
These damages are mainly formed by DSB (Double- Strand Break) lesions 
[9]. Around these DSBs it will be formed repair complexes. These complexes 
include at eukaryotes the proteins: Mrell, Rad50 and NBSl (MRN complex) 
and they are the signal for activation of ATM. The DSBs repair protein 
complexes count only for the initial activation of ATM, because the main 
activation is an autocatalytical process. ATM, at its own, represents the 
signal for activation of gene p53. Depending on the level of ATM, p53 will 
lead to two outcomes for the cell: one is the cell cycle arrest induced by a 
low level or a brief elevation of p53 protein, and the other is the apoptosis 
induced by a high level or a prolonged elevation of p53 protein [8] . Each of 
these two outcomes could not be an option for all the cells. For example, 
apoptosis is not accepted for neurons or myocardial muskular cells because 
they do not divide in adult life, so these cells will choose for cell cycle arrest. 
On the other hand for the enterocytes (cells of digestive tube) apoptosis is a 
common option, because these cells divide themselves very quickly and their 
lifetime is no longer than 2 days. Now is clear that should be a very good 
control of p53 activity in such a manner that the cell goes on right pathway 
(i.e. apoptosis or cell cycle arrest). This control is achieved with the help of 
mdm2 gene with which p53 makes a feedback loop [8, 12]. 

In the last years, the approaches of P53 dynamics as response to DNA 
damage comprise modelings in which are described three distinct subsys- 
tems: a DNA damage repair module, an ataxia telengiectasia mutated (ATM) 
switch and the P53-Mdm2 oscillator. 

In what follows we will consider a model only for the third module. The 
variables of the model are: Xi P53-mRNA concentration, X2 Mdm2-mRNA 
concentration, yi P53-protein concentration and ya Mdm2-protein concen- 
tration. 

We consider P53-Mdm2 model with kernel delay given by: 



2 



Xi{t) = ai - a2Xi{t), 

yi{t) = bixi{t) - b2yi{t) - 612^1 (t)y2(t), ^-^^ 

i2{t) = /o°° h{s)f{yi{t - S))ds - C2X2{t), 

M) = /o°° k2{s)x2{t - s)ds - d2y2{t) - di2yi{t)y2{t) 

where: 02, C2 are the rates for mRNA degradation, 62, c^2, ^12, du are the rates 
for proteins degradation. The function / : ]R+ — > ]R+, is the Hill function, 
given by: 

f{x) = ^ 

with e INT*, a > 0. The parameters ai, 02, 61, 62, C2, 612, (^2 <^i2 of the 
model are assumed to be positive numbers less or equal to one, the memory 
functions ki, k2 that reflect the influence of the past states on the current 
dynamics are a nonnegative bounded function defined on [0, 00) and 

/ ki{s)ds — 1, / ski{s)ds < 00, i — 1,2. 
Jo Jo 

The memory function is called delay kernel. The delay becomes a discrete 
one when the delay kernel is a delta function at a certain time. Usually, we 
employ the following form: 

kiis)^^sPe-''*', qi>0,i = l,2,p>0 

for the memory function. When p — and p — 1 the memory function are 
called "weak" and "strong" kernel respectively. 

For ki{s) — 5{s — Ti), k2{s) — S{s — T2), ti > 0, T2 > the model is given 

by: 

xi(t) = ai - a2Xi(t), 

yi{t) = bixi{t) - b2yi{t) - 6i2?/i(i)?/2(i), 

^2(0 = f{yi{t - ^1)) - C2X2(t), 

Mi) ^X2{t- T2) - d2y2{t) - di2yi{t)y2{t). 

In (2) with Ti = T, T2 = 0, di2 = we obtain the model from [11] and it 
was studied in [10] which suggests that there is an oscillatory behavior based 
on observations obtained using only numerical simulations. 
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In present paper we will analyze the model (1) with di2 = with the 
following initial values: 

X2{0) = MO),ee (-00,0], 1/2(0) = y2, 

with xi > 0, ^2 > 0, (pi{9) > 0, (^2(9) > 0, for all 9 e (—00, 0] and (p2 are 
different iable functions. 

Also, two delays appear in the leukopoiesis model which is analyzed in 

[!]• 

The paper is organized as follows. In section 2, we discuss the local 
stabihty for the equihbrium state of system (1), with di2 — 0, if the delay 
kernels ki, k2 are delta functions, and k\ is delta function, k2 is weak kernel 
function. We investigate the existence of the Hopf bifurcation with respect to 
the parameters of the delay kernels fci, k2. In section 3, the direction of Hopf 
bifurcation is analyzed by the normal form theory and the center manifold 
theorem introduced by Hassard [4] . Numerical simulations for justifying the 
theoretical results are illustrated in section 4. Finally, some conclusions are 
made. 



2. Local stability and the existence of the Hopf bifurcation. 

We consider the model: 

Xi{t) = tti — 02X1 (t), 

ijiit) = bixi{t) - b2yi{t) - hi2yi{t)y2{t), 
X2{t) = /o°° ki{s)f{yi{t - s))ds - C2X2(t), 
y2{t) = /o°° k2{s)x2{t - s)ds - d2y2{t). 

Proposition 1. If < bi and yio e (0, ^* solution of equation 

a2t>2 

- /3a;" + 7a; - 5 = 

where 

a = 02(612 + b2C2d2),P = aibiC2d2, 7 = a2&2C2C?2a", S = Oi&i 02^20" 
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then the equilibrium point X* of system (3) has the coordinates: 

01 aihi - 0262^10 

Xio = — , X2Q = rf2?/20, ?/20 = T ■ 

02 Oi2yioa2 

We consider the following translation: 

= Ml + xio, yi = 1*2 + 2/10, a:;2 = M3 + 2:20, 2/2 = 1*4 + y2o- (4) 
With respect to (4), system (3) can be expressed as: 
ui{t) = -02^1 (t), 

U2{t) = - (62 + bi2y20)u2{t) - bi2yioU4{t) - bi2U2{t)u4{t) , 

Mt) = So" ki{s)f(u2(t - s) + yio)ds - C2{u3{t) + X20), 

U4{t) — /q°° k2{s)u3{t — s)ds — d2U4{t). 

System (5) has (0, 0, 0, 0) as equihbrium point. 

To investigate the local stability of the equilibrium state we linearize 
system (5). We expand it in a Taylor series around the origin and neglect 
the terms of higher order than the first order for the functions from the right 
side of (5). We obtain: 

U{t) = AU{t) + BiUi{t) + B2U2{t), (6) 

where 
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with p = /'(yio), U{t) = {ui{t)Mt)Mt)Mt)f, U,{t) = Pi{s)u,{t - 
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s)ds,J^ Pi{s)u2{t - s)ds,j^ pi{s)uz{t - s)ds, Pi{s)u4{t - s)ds)'^, i= 1,2. 

The characteristic equation corresponding to system (6) is A(A) = 0, 
where 

A{X)^det{XI-A-{ ki{s)e-^'ds)Bi- { k2{s)e-^'ds)B2). (9) 

Jo Jo 

From (7), (8) and (9) it results: 

A(A) = (A + a2)Ai(A) 

where 

Ai(A) = A=' + j92A'+PiA + po + r(/ h(s)e-^'ds)( k2(s)e-^'ds) (10) 

Jo Jo 

with 

P2 = f'2 + C2 + C?2 + fcl2l/20,Pl = (C2 + C?2)(&2 + &122/2o) + C2C/2 /-^-^n 

Po = 02^2(62 + bi2y2o),r = pbuVw 

The equilibrium point X* — {xw, yw, X20, 1/20)^ is locally asymptotically 
stable if and only if all eigenvalues of A (A) = have negative real parts. 
Because 02 > 0, we will analyze the equation Ai(A) = 0. The analysis of the 
sign of real parts of eigenvalues is complicated and a direct approach cannot 
be considered. 

We will analyze the eigenvalues for the equation Ai(A) = if the delay 
kernels ki and k2 are delta functions or ki is delta function and k2 is weak 
function. 

Proposition 2. // ki{s) = S{s — ti), k2{s) = 6{s — T2), n > 0, ri > 
then: 

(i) function (10) is given by: 

Ai(A, r) = A^ + P2A' + piX + P0 + re-^^ (12) 

where t — ti + T2. 

(ii) if Ti = 0, T2 = then the equilibrium state X* of system (5) is locally 
asymptotically stable if and only if 

PiP2>Po + r (13) 
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where Pi, p2, Po, r are given by (11). 

We are looking for the values tq so that the equilibrium point X* changes 
from local asymptotic stability to instability or vice versa. If the relation 
(13) holds, for T > sufficiently small all roots of Ai(A,r) = have nega- 
tive real parts. The critical delay Tq is the smallest positive value of r where 
Ai(A,r) = has imaginary roots. Let A = ±za; be these solutions with 
a; > 0. Separating real and imaginary parts of Ai(za;, r) = we obtain: 

rcos{ujT) — p2ijJ^ — Po, rsin{uT) — ujp\ — uj^ . (14) 

A solution of (14) is a pair (a;o, To)where uq is a positive root of the 
equation: 

+ {pI - 2pi)x^ + {pj - 2pqP2)x^ + - = 
and To is given by: 

Tq — Uii Ciy 7^ . 

^0 -P2'^0 + Po 

Prom (12), we obtain: 

y _ (iA _ Ar 

~ rfr ~ ~e^^(3A2 + 2p2A + pi) -rr' ^ ' 

Then, we evaluate (15) at A = iujq and r = Tq and obtain: 
where 

h = {Pi - Suq)cos{ujoTo) - 2p2UJosin{uJoTo) - rro 
k = {Pi - Suj^)sin{ujoTo) + 2p2^ocos(a;oro). 

From the above analysis and the standard Hopf bifurcation theory [4], we 
have the following result: 

Proposition 3. If pi, P2, Po, f satisfy (13) and po < r, for r — To, 
CO — u>o then: 

and a Hopf bifurcation occurs at the equilibrium state X* as r passes through 
To- 
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Proposition 4. If ki{s) = S{s — ti), k2{s) = q2e ^^^^ n > 0, ^2 > then: 

(i) function (10) is given by: 

Ai (A, n) = (A + q2) (A^ + p2X^ + pi A + po) + rgse'^^^ ; (16) 

(ii) ifTi — 0, then the equilibrium state X* of system (5) is locally asymp- 
totically stable if and only if 

D2 = (P2 + q2){pi + q2P2) - {po + Piq2) > o, 

D3 = {po + Piq2)D2 - {P2 + q2f{q2PQ + rq2) > 0. 



We are looking for the values t^q so that the equihbrium point X* changes 
from local asymptotic stability to instability or vice versa. The critical delay 
r*Q is the smallest positive value of ri where Ai(A,rj*Q) = has imaginary 
roots. Let A = ±2cj be these solutions with c<; > 0. Separating real and 
imaginary parts of Ai(ia;, r) = we obtain: 



q2rcos{ujTi) = -uj"^ + {pi + ^2^2)^^ - ?2Po, 
q2rsin{uTi) = -{p2 + ^2)^^^ + {Po+Piq2)^- 



(18) 



A solution of (18) is a pair (cjio, rj*Q)where uj^q is a positive root of the 
equation: 



where 



+ nix^ + n2X^ + uzx^ + 77.4 = 



ni = (p2 + q2f - 2(pi + ^2^2), 

n2 = {Pi + q2P2f + 2g2Po - 2(po + Piq2){P2 + ^2), 

ns = iPo+Piq2f - 2g2Po(Pi + g2^'2), 

22 22 

n4 = Pt,q2-r qi- 



and T^Q is given by 

T*o = —arctg 

Prom (16), we obtain: 
dX 



X 



{P2 + q2)cOw - (Po + Plg2)'^10 

^10 - (Pl + q2P2)^W + ?2P0 ' 



Arg2 



(19) 



dri e^^i (4A3 + 3(p2 + q2)X'^ + 2(pi + q2P2)X + po+ ^1^2) - rq2Ti 



(20) 
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Then, we evaluate at A = iojio and ri = t^q and obtain: 

x'/ * ^ _ cc;iorg2^2o . , u)iorq2ho 
no "T '20 no "T '20 

where 

^10 = (-3(P2 + q2)^lo +P0+ Piq2)cOs{uJiQTiQ) 

+(4a;?o - 2(pi + q2P2)^iQ)sin{LOwTw) - rr-^qi 

ho = (-3(j02 + g2)ti^?o + + Pig2)s^?T-('^loTlo) 

+(4a;^o + 2(pi + ?2P2)ti^io)cos(a;oro). 
Wc have the following result: 

Proposition 5. If pi, p2, po, r, 52 satisfy (17) for ti = t^q, lv = ujiq 
then: 

V""^! / A=ja;io,Ti=r*o 'lO + '20 

and a Hopf bifurcation occurs at the equilibrium state X* as Ti passes through 

T* 

no- 



3. Direction and stability of the Hopf bifurcation 

In what follows, we will study the direction and stability in two cases: in 
the first case the both kernels are delta function and in the second case the 
kernel ki is delta function and the kernel k2 is weak function. 

3.1. The case ki{s) = 5{s — ti), k2{s) = 5{s — T2), ti > 0, T2 > 0. 

For ki{s) — S{s — Ti), k2{s) — 5{s — T2), ti > 0, ti > from 
Proposition 2, wc obtained some conditions which guarantee that system 
(5) undergoes Hopf bifurcation at r = Tq. In this section, we study the 
direction, the stability and the period of the bifurcating periodic solutions. 
The used method is based on the normal form theory and the center manifold 
theorem introduced by Hassard [4] . We know that if r = Tq then all roots of 
^i(^) ''0) = 0, when Ai(A, Tq) is given by (12), other than ±iujQ have negative 
real parts and any roots of the form A(t) — q:(t) ± ia;(r) satisfies Q;(ro) = 0, 

"^("^o) — <^o ^-nd ^'^1'^°'^ 0. 

dr 

Suppose that for given ai, 02, hi, 62, &12, C2, a, c?2 there is tq for which 
Ai(A, To) = exhibits a Hopf bifurcation. We consider rio = tq — T2, where 
T2 < 2to and Ti = Tiq + G H. We regard n as the bifurcation parameter. 
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For $ e = C([— Ti, 0],^^"^) we define a linear operator: 
L^($) = A$(0) + 5i$(-ri) + S2$(-r2) 
where A,i?i, i?2 are given by (7), (8) and a nonlinear operator: 

F{pi, $) = (0, -6i2$2(0)$i(0), ^/92$^(-Ti) + ^p,^l{-r^),Of + 0(|$|^) 

where $ = ($i, $2, ^3, P2 = /"(yio), Ps = /"(z/io)- 

By the Riesz representation theorem, there exists a matrix whose com- 
ponents are bounded variation functions r]{0,n) with 9 e [— Tio,0] such that 

r d?7(^,/i)$(^), [-rio,0]. 
-'—no 



A, 9^0 

B26{e + T2), 9e[-r2,0) 
B^5{9 + t{), [-rio,-T2). 




0, ee[-rio,0) 

Then, we can rewrite (5) in the following vector form 

iit = A{ij)ut + R{y)ut (21) 

where Ut — u{t + 9), iox 9 & [— tio, 0]. 

As in [4] the bifurcating periodic solutions u{t, f^) of (21) are indexed by 
a small parameters e, e > 0. The solution u{t,ii{e)) has amplitude 0{e), 
period T{e) and nonzero Floquet exponent (3{e) with f3{0) — 0, where under 
our conditions /i, T and /3 have convergent expansions: 



10 



For e C^([0,rio],(D**^), the adjoint operator ^* of A is defined as: 



A*^{s) = { ds 



s e (O,rio] 



For $ e C([-Tio,0],(D^) and * e Ci([0, Tio],(D*^) we define the foUowing 
bihnear form: 

< $(^) >= ^(o)^$(o) - r r _ (22) 

where 1]{9) — r]{9, 0). 

Then, it can verified that A* and A are adjoint operators with respect to 
this bihnear form. 

For system (21) we have: 

Proposition 6. // Ai = icoo, A2 = Ai then: 

(i)The eigenvector of A{0) corresponding to Ai is 

h{e)=ve^'\ ^e[-Tio,0] 
where v = (^^1,^2, v^, v^Y , 

vi = 0,^;2 = -(Ai + d2)(Ai + 02), V3 = -pe^'^'°{\i + d2),Vi = -pe^'''\ 

To = Tio + T20- 

(a) The eigenvector of A* corresponding to A2 is 
h*(s)^we^'', se[0,oo) 
where w = {wi,W2, w^, w^)'^, 

a2 + X2 6^1^2^12^10(02 + A2) 

Ol (C2 + A2)(a2 + A2jOl 

fci2t/io(a2 + A2)e^i"2 

^4 = — r"71 TTT TT^ 

0i(a!2 + A2)(c2 + d2) 
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Ol (C2 + A2)(d2 + A2)6i 

(~ Air,- ^ ^12^10 (02 + A2) 

01(02 + A2j 

(in) With respect to (22) we have: 

<h*,h>=l, < h*,h>=<h*,h>=0, <h*,h>=l. 

Using the approach in [2], we next compute the coordinates to describe 
the center manifold fio at /x = 0. Let Ut = u{t + 0),9 G [— rio,0), be the 
solution of system (21) when /j, — 0. 

We define 

z(t) =< h\ Ut >, w{t, 0) = ut{e) - 2Re{z{t)h{9)). 
On the center manifold flo, we have: 

w{t,e) ^w{z{t),z{t),9) 

where 

w{z, z, 9) = ^20(6")— + wii{9)zz + ^02(6')— + ^30(6*)-^ + ■ • • 
2 2d 

in which z and z are local coordinates for the center manifold VLq in the 
direction of h* and h* and wq2{.9) = W2o{9). 

For solution Uf G of equation (21), as long as ^ = 0, we have: 

z{t) = Xiz{t) + h*{0)F{w{z{t),z{t),Q) + 2Re{z{t)h{Q))) = 
Xiz{t) +g{z,z) 



where 



z^ z^ z^z 



9{Z: Z) = g20— + guZZ + ^02 Y + ^21— + • • • 



Proposition 7. For the system (21) we have: 

g20 = -2bi2V2ViW2 + p2^^iw3e2^2rio^ 

gn = -bl2iv2V4 + V2Vi)w2 + P2V2V2'W3, (23) 

go2 = -2bi2V2ViW2 + p2vlw2,e^^'-'^^° , 
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(u) 

/Q\ Aie 9ii Xoe I 771 

wu{0) = —ve^ — —ve^ + E2, 
Ai Ai 

where Ei = {Eu, ^21, ^31, £^41)^ and E2 = (-E12, E22, E32, E42) 

P P 

E31 = (2Ai + rf2)e-2A2r2^^^ 

P2wi(2Ai + 62 + 6121/20) - 26i2f2f4P 

-c/41 — 



P62yi0 + (2Ai + 62 + ^'l2?/20)(2Ai + C2)e-2^2ro 

P HP ^2-t^2^2 C20?2 4 4 
-C/12 = U, £/22 = \ -C/2 5 -^32 = "2-C/2 

P&12(''2''1 + ''2''l) - P2l'2l'2Q>2 + h2lj2o) 
-C/42 = 



(^'2 + 6l22/20)c2C?2 + Pb2yw- 



(in) 



§21 = -3612(^210420(0) + 2t;2ty41l(0) + 1^4^220(0) + 2W21l(0)f4)w2 + 

+W3[6p2{2v2e^^^^ - W2iii-Ti) + 6i;2e^i^i«;22o(-ri)) + 3p3vie^^^^^V2e^^^^], 

(24) 

withw2o{0)={wi2o{0), W22o{0), W32o{0), W42o{0)V and Wn{9) ^{win{9) , W2n{0), 

Based on the above analysis and calculation, we can see that each gij 
in (23), (24) are determined by the parameters and delay from system (21). 
Thus, we can explicitly compute the following quantities: 

^^i(0) = ^^(920911 - 2\gn\' - ^1^02!') + ^ 

ReTdm ^ Jm(Ci(0))+//2/my(0) . _ 
= - Rem ' = ^ 2i?e(Ci(0)), 

(25) 

where A'(0) is given by 



A'(0) 



^e^^(3A2 + 2p2A+pi -r) 



X=iU!o,T=To 



13 



In summary, this leads to the following result: 

Theorem 1. In formulas (25), fi2 determines the direction of the Hopf 
bifurcation: if H2 > 0(< 0), then the Hopf bifurcation is supercritical (sub- 
critical) and the bifurcating periodic solutions exist for r > ro(< Tq); P2 
determines the stability of the bifurcating periodic solutions: the solutions 
are orbitally stable (unstable) if P2 < 0(> 0); and T2 determines the pe- 
riod of the bifurcating periodic solutions: the period increases (decreases) if 
T2 > 0(< 0). 

3.2. The case ki{s) = 6{s - n), k2{s) = q2e-'^^% n >0, q2> 0. 

For ki{s) = 6{s — ri), k2{s) = §26"^^*, ti > 0, §2 > 0, system (5) is 
given by: 

ui{t) = -a2Ui(t), 

il^{t) = biUi{t) - (62 + hi2y20)u2{t) - bi2yi0Uiit) - bi2U2{t)Ui{t), 

ilsit) = f{u2{t - Ti) + ?/io) - C2{u3{t) + X20), (26) 

U4{t) = U5{t) - d2U4{t), 

U5{t) = q2{u3{t) - U5{t)). 

We expand it in a Taylor series around the origin and neglect the terms 
of higher order than the first order for the functions from the right side of 
(26). We obtain: 

U{t) = A,2U{t) + B,2U{t-n), 

where 



A 



12- 



/-a2 








\ 


bi -(62+6122/20) 





-612I/10 








-C2 
















d2 


1 


^ 







-Q2 J 




/ 





^ 





















B12 — 


p 





























^ 





°J 







(27) 



(28) 
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with U(t) = {ui(t),U2(t),U3(t),U4(t),U5{t)), U{t-Ti) = {ui{t - Ti) , U2{t - 
Tl),U3{t - Ti),U4{t - Ti),U5{t - Ti)Y . 

Let tIq given by (19) and Ti = t*q + e H. We regard ji as the 
bifurcation parameter. 

For $ e = C^([— Ti, 0],(E^) we define a hnear operator: 



where $ = ($i, $2, $3, *4, ^s)"^, p2 = /"(yio), Ps = /"(yio)- 

By the Riesz representation theorem, there exists a matrix whose com- 
ponents are bounded variation functions r]{6,n) with 9 e [— tio,0] such that 



^i2^(<f ) = ^12$ (0) + 5i2<f (-ri) 



where A12 ,Bi2 are given by (27), (28) and a nonhnear operator: 



Fi2(//, $) = (0, -6i2$2(0)$i(0), -p2$^(-ri) + -p3$^(-ri), 0, 0)^ + 0(|«|^) 




We can choose 



A 



12, 



^ = 



Bi2S{9 + Ti), ^e[-rro,0). 



For $ e C we define: 





Then, we can rewrite (26) in the foUowing vector form 



Ut = Al2{lJ')Ut + Ri2{fJ')Ut 



(29) 



where Uf = u{t + 9), for 9 e [-t^q, 0]. 
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For e C^([0,rio],(D*^), the adjoint operator AI2 of ^ is defined as: 

[ J^,,Jv''{t,0)^{-t), s = 0. 

For $ e C([-Ti*o,0],(D^) and * e Ci([0, Ti*o],(D*^) we define the following 
bilinear form: 



< ^9) >= ^{Ofm - f f ^^(C - e)drii2{e)^i)di, (30) 



where r7(^) = 'q{9, 0). 

Then, it can verified that a-nd are adjoint operators with respect 
to this bilinear form. 

For system (29) we have: 

Proposition 8. // Ai = iojiq, A2 = Ai then: 

(i)The eigenvector 0/^12(0) corresponding to Ai is 

h{e)=ve^^\ ^e[-r*o,0] 

where v = {vi,v-2, v^, V4, v^Y , 

r, (Ai+g2)(Ai+C2) q2 

Vi = Q,V2 = e 1 = Ai + g2,V4 = , , , > '^5 = ^2- 

Pi Ai+d2 

(a) The eigenvector of A\2 corresponding to A2 is 
h*{s)^we^'% se[0,oo) 
where w = {wi, W2, W3, W4^, w^y , 

h 1 



(A2 + 02)77'^^ ?7'^^ (c2 + A2)(rf2 + A2)(g2 + A2)?7' 

_ ^12^/10 012^10 

(c?2 + A2)?7' (c?2 + A2)(g2 + A2)?7 



{C2 + X2)id2+X2){q2 + X2f' {d2 + X2)'" (^^2 + A2) (^2 + A2) 



(c2 + A2)(rf2 + A2)(Q2 + A2)Ai 
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(in) With respect to (20) we have: 

<h*,h>^l, <h*,h>^<h*,h>^0, <h*,h>^l. 



Using the approach in [2], we next compute the coordinates to describe 
the center manifold Qq at n = 0. Let Ut — u{t + 0),9 e [—Tio,0), be the 
solution of system (29) when /i = 0. 

We define 

z{t) =< h*, Ut >, w{t, e) = ut{e) - 2Re{z{t)h{e)). 
On the center manifold Q,o, we have: 

w{t,0) = w{z{t),z{t),0) 

where 

^2 Z^ 

W{Z,Z,6) = W2o{9)— + Wii{9)zZ + Wo2{0)— + W3o(6')— + . . . 

in which z and z are local coordinates for the center manifold Qq in the 
direction of h* and h* and Wo2{d) = 'W2o{d)- 

For solution ut G flo of equation (29), as long as = 0, we have: 

z(t) = Xiz{t) + h*{0)F(w(z(t),z(t),0) + 2Re(z{t)h(0))) = 
Xiz{t) + g{z,z) 



where 



z^ z^ z^ z 

giz, z) = g2o— + guzz + go2— + ^21— + • • • 



Proposition 9. For the system (29) we have: 

g2o = -2bi2V2ViW2 + p2vlw^e^^''''*o ^ 

gu = -bl2{v2V4 + V2Vi)w2 + P2V2V2W3, (31) 

go2 = -2bi2V2V4W2 + p2vlw2,e^^^'^w ^ 



(u) 



920 -x,e 6*02 - , 771 2Aie 
W2o{0) — — —ve ^ — —T-ve^ + Eie^ 

Ai _ 3Ai 

wii{9) = —ve'^^" — —ve''^" + E2, 
Ai Ai 
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where Ei = {Eu, E21, E31, E41, E^if and E2 = {E12, E22, ^32, E42, E^2f 

jj, ^ 0^22-^220 ~ «12-^320 j-^ ?2 + 2Ai 

= U, £/2l = ; , -C/31 — -C/51, 

011022 + 012021 q2 

J-, _ 1 771 171 O21-F22O + 0,11^320 

-C/41 — — — --— £^51, -t/51 — : 

O2 + 2X1 011022 + 012021 

On = 2Ai + 62 + O12Z/10, ai2 = , , ^, , 

02 + 2Ai 

2^^^. (c2 + 2Ai)(g2 + 2Ai) 

Q2 

F220 = -2612 W2W4, F320 = p2^^|e2^2r*o^ 

p - n p - "^22^^211 - C12F3II p _ p p _ 1 p 

£/12 — U, i^22 — ; 5 -C'32 — -C/52, -C'42 — ^-£^52, 

CiiC;i2 + C12C21 (22 
C2l-r211 + C11F311 

-£^52 — ; 

C11C22 + C12C21 

cii = 62 + 6i2?/2o, C12 = buy 10, C21 = pi, C22 = C2 

i^211 = -&12(^^2^^4 + ^^2^^4), i^311 = ^^2^^2/02- 



(in) 

921 = -36l2(^^2«'42o(0) + 27J2W'41l(0) + t^4«'22o(0) + 2«;211 (0)7J4) W2 + 

^3 [6p2 {2v2e^'^^o - W2n (-rj^o) + 6^26^1^10 W220 (-^1*0) ) + 3p3vy^^^^ov2e^^^^, 

(32) 

withw2o{9) = {wi2o{0), W22o{0), W32q{0), W42o{0), Wb2o{0)) and wii{9) ={wiii{6) , 

W2ii{e), Willie), Willie), w^nie)), e e [-ti*o,o]. 

Based on the above analysis and calculation, we can see that each gij 
in (31), (32) are determined by the parameters and delay from system (26). 
Thus, we can explicitly compute the following quantities: 

Ci(0) = ^^{920911 - 2\gii\' - ^\go2\') + ^ 

RefCim ^ /m(Ci(0))+p2/mA-(0) ^r,.r((^^^ 
^ - Rem '^^ ^ Wo ^ 2i?e(Ci(0)), 

(33) 

where A'(0) is given by 
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eAri(4A3+3(p2+g2)A2+2(pi+g2P2)A+po+Pi92)-rg2'ri 



Arg2 



>F=ja;io,Ti=7-*Q 



We have: 



Theorem 2. In formulas (33), 112 determines the direction of the Hopf 
bifurcation: if 112 > 0(< 0), then the Hopf bifurcation is supercritical (sub- 
critical) and the bifurcating periodic solutions exist for t\ > t*o(< Tiq); 

determines the stability of the bifurcating periodic solutions: the solutions are 
orbitally stable (unstable) if ^2 < 0(> 0); andT2 determines the period of the 
bifurcating periodic solutions: the period increases (decreases) ifT2 > 0(< 0). 

4. Numerical examples. 

For the numerical simulations we use Maple 9.5. In this section, we 
consider system (6) with ai = 2,02 = 0.55, hi = 1, 62 = 0.8, C2 = 0.1, 
612 = 1.5, d2 = 0.1, a = A, n — 2. We obtain: xio — 3.636363636, yio = 
0.8347719895, ?/2o = 2.370744013, 0:20 = 0.2370744013. 

In the first case, ki{s) — 5{s — ri), ^2(5) = (^(s — 7-2), for r2 = 3, we 
have: Wio = 0.1324013896, 112 = -0.4204703301, /32 = 0.2799153884, T2 = 
0.0005051758260, To = 9.541873607 Then the Hopf bifurcation is subcritical 
and the bifurcating periodic solutions exist for r > Tq; the solutions are 
orbitally unstable and the period of the solution is increasing. The waveforms 
are displayed in Figl and Fig2 and the phase plane diagrams of the state 
variables yi{t), 2/2 (^) and yi{t — r), yi{t) are displayed in Fig3 and Fig4: 
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Fig.3. {yi{t),y2{t)) 



Fig.4. (yi(t-T),yi(t)) 





1.834 1.8344 1.E 



1.834 1.8344 1.8348 1.8352 



In the second case, ki{s) = S{s — ti), k2{s) 
have: = 0.1290621026, ^2 = 



q2e' 



'q2s 



for q2 = 0.5, we 



-0.5993860816, /?2 = -0.7476750590, T2 = 
0.1798944390, t*q = 32.37014890. Then the Hopf bifurcation is subcritical 
and the bifurcating periodic solutions exist for ti > t^q] the solutions are 
orbitally stable and the period of the solution is increasing. The waveforms 
are displayed in Fig5 and Fig6 and the phase plane diagrams of the state 
variables yi(t), y2(t) and yi{t — r), yi{t) are displayed in Fig7 and Fig8: 
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Fig.5. {t,yi{t)) 



Fig.6. {t,y2{t)) 




5. Conclusions. 



As in our previous models [11,13], we obtain an oscillatory behavior simi- 
lar to that observed experimentally [5]. The conclusion is not surprising, but 
is useful as this model provides a more accurate approach of the interaction 
P53-Mdm2. 
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The improvements of the model from [11] done in the present paper, 
proved usefulness as we obtained a smoother modelling of the phenomenon 
and the oscillating behavior remained which as similar with that from [5] . 

Using the method from this paper, we will do a qualitative analysis of 
the model from [9] in our future papers. 
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